Entanglement swapping with local certification: 
Application to remote micromechanical resonators 
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We propose a protocol for entanglement swapping which involves tripartite systems. The genera- 
tion of remote entanglement induced by the Bell measurement can be easily certified by additional 
local measurements. We illustrate the protocol in the case of continuous variable systems where 
the certification is effective for an appropriate class of three-mode Gaussian states. We then apply 
the protocol to optomechanical systems, showing how mechanical entanglement between two remote 
micromechanical resonators can be generated and certified via local optical measurements. 
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The protocol of entanglement swapping fll| is recog- 
nized as a fundamental tool for quantum communication. 
It is one of the central mechanisms at the basis of quan- 
tum repeaters ,2;] , which allow to distribute entanglement 
at long distances in an efficient manner. Very recently, 
the importance of entanglement swapping has been also 
recognized in quantum cryptography, where this tech- 
nique introduces a Hilbert space filtering able to exclude 
side-channel attacks 0). 

Entanglement swapping allows to entangle two distant 
systems which never interacted. The standard protocol 
considers two distant users, Alice and Bob, each pos- 
sessing a bipartite entangled system, RiBi for Alice and 
R2B2 for Bob. The users send their systems Bi and B2 
to an intermediate station (Charlie). Here an appropri- 
ate Bell measurement is made over the B systems and 
the outcome communicated back to Alice and Bob. As 
a result, the two remote systems i?i and R2 become en- 
tangled. 

After the completion of this procedure, an important 
point is the experimental verification of the generated 
entanglement, which generally involves direct measure- 
ments on the two remote systems, i?i and i?2j followed by 
classical communications. This is fine as long as systems 
i?i and R2 are easily measurable, which is not always the 
case though. For instance, an important scenario is repre- 
sented by the entanglement swapping of optomechanical 
systems, where i?i and R2 are mechanical modes while 
Bi and B2 are optical modes In this case, entan- 
glement generation is easy, since the Bell measurement is 
done on the optical modes, but entanglement verification 
is much more difficult, since it involves quantum- limited 
measurements on mechanical modes [5|-l9|, furthermore 
placed at distant nodes. 

Whenever the verification of entanglement involves dif- 
ficult quantum measurements, an experimentalist may 
therefore consider an alternative approach which is based 
on the local certification protocol proposed in this Letter. 
The key idea is to extend the initial systems of Alice and 
Bob from bipartite to tripartite. We add a certification 



system (C) to each site, so that we have tripartite sys- 
tem RiBiCi for Alice and another R2B2C2 for Bob. In 
this new protocol, also the certification systems Ci and 
C2 are sent to Charlie's station, where they are locally 
measured. This additional measurement enables Char- 
lie to certify the generation of entanglement for the two 
remote systems Ri and R2. 

In practical implementations, the certification systems 
are chosen to be easy to detect (e.g., ancillary optical 
modes in the case of optomechanical systems). However, 
in order to achieve certification, we cannot use arbitrary 
tripartite states but a suitable class of states must be 
engineered, that we call "certifying" . In our derivation, 
where we consider continuous variable systems, we show 
that these states correspond to a generic class of three- 
mode Gaussian states satisfying specific conditions of pu- 
rity. 

As an application of our protocol, we consider the 
scenario where Alice and Bob possess two identical op- 
tomechanical systems, each formed by two output modes 
of an optical cavity coupled to a micro-mechanical res- 
onator. One traveling optical output mode from each 
site is mixed and employed for the Bell measurement 
generating entanglement between the two remote micro- 
mechanical resonators, while the ancillary optical output 
mode from each site is used for entanglement certifica- 
tion. The protocol proposed here is easier to implement 
than a standard entanglement swapping protocol based 
on bipartite optomechanical states at Alice and Bob sites. 
In fact, the preparation of the two initial certifying states 
is not more difficult than preparing bipartite optome- 
chanical entanglement, because adding a second driven 
cavity mode is a standard option in current cavity op- 
tomechanics ?|. Instead, the experimental verification 
stage is much simpler because it involves only local opti- 
cal measurements at Charlie's site, avoiding any nonlo- 
cal measurement on the mechanical resonators. We show 
that using state-of-the-art apparatuses, it is possible to 
generate robust entanglement between two remote meso- 
scopic systems, which can be easily certified by local opti- 
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cal measurements. Such a scheme could be exploited for 
fundamental tests of quantum mechanics at the macro- 
scopic level U I ■ 

The Letter is organized as follows. We start with the 
definition of tripartite certifying states. Then, we de- 
scribe the entanglement swapping protocol with local cer- 
tification. Finally, we apply the protocol to the case of 
optomechanical systems. 

Tripartite certifying states. Consider a bosonic sys- 
tem which is composed of three modes, labelled as R 
(for remote mode), B (for Bell mode), and C (for cer- 
tification mode). Each mode is described by quadra- 
ture operators Xk and pk, with commutation relations 
[xk,Pk'] = iSkk' and [xk,Xk>] = \pk,Pk'] = 0. We as- 
sume that the system is prepared in a zero-mean Gaus- 
sian state. This means that its characteristic function is 
Gaussian $(k) = exp(-k"rVk/2), where k e r^.^^ y 
is the covariance matrix (CM). We can always write this 
CM in the form 



V = 



RDF 

DT B E 
pT C 



(1) 



where the blocks R, B, C, D, E, F are 2x2 submatrices. 
By assuming non-singular blocks, this CM can always 
be simplified to the standard form of Ref. via local 
unitary operations. We implicitly assume this standard 
form, where R = rl, B = 61, C = cl (I is the identity 
matrix) and D, E are diagonal [l3|. From the tripartite 
CM, we can extract the two submatrices 
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where Vrb describes the remote-Bell modes, and Vbc 
the Bell-certification modes. Denoting with /i^g = 

[4VdetVRB]"\ /^ec = [4VdetVBc] 



and 



Mb 



2VdetB the purities of the RB, BC, and B sub- 
systems respectively, we say that the 3-mode Gaussian 
state is certifying when 



Mrb > A^BC > Mb- 



(3) 



Notice that the certifying condition implies that both 
RB and BC subsystems are entangled. In fact, 

whenever 
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Mbc ^ Jin, 
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1. Therefore, 
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for a certifying state one has 



/i^^^ which, using the result of 



15| . implies that the BC subsystem is entan- 
gled. The same argument applies to the RB subsystem. 

Entanglement swapping with certification. Consider 
two distant users, Alice and Bob. Each user has a sys- 
tem composed of three bosonic modes: these are R^BiCi 
for Alice and R2B2C2 for Bob (see Fig. [T^). For sim- 
plicity we consider here the case where the two sys- 
tems are prepared in the same certifying state with CM 





FIG. 1: (Color online) (a) Entanglement swapping with, local 
certification. Alice and Bob possess two tripartite systems, 
RiBiCi and R2B2C2, prepared in a certifying state. Alice 
and Bob keep their remote modes Ri and R2, while sending 
the other modes to Charlie. Bell modes Bi and B2 are sub- 
ject to a Bell measurement with outcomes {x-,p+} stored 
by Charlie. Certification modes Ci and C2 are then subject 
to local measurements. The detection of C1C2 entanglement 
certifies that potential entanglement is also present between 
between Ri and R2. This remote entanglement is then estab- 
lished by the communication of {x-,p+} to Alice and Bob. 
(b) Optomechanical system mounted on each site. The remote 
mode is the mechanical resonator (quadratures q and p) while 
the Bell and certification modes are the two output optical 
modes with annihilation operators 5^°' and a^"', respectively. 
These modes derive from corresponding intracavity modes at 
and Oc coupled with the mechanical resonator. 



VrjBiCi = VR2B2C2 ~ "^i which can be expressed in the 
blockform of Eq. ([T|) . 

In the first step of the protocol, Alice and Bob retain 
their remote modes while sending Bell and certification 
modes to Charlie. In other words, Alice keeps mode 
Bob keeps mode R2, while Charlie gets modes Bi, B2, 
Cl and C2 (see Fig.[T}i). 

In the second step, Charlie performs a CV Bell mea- 
surement on the Bell modes Bi and B2 by using a 
balanced beam-splitter and two homodyne detectors. 
This corresponds to measuring the two combinations of 
quadratures = XB2 ~ ^Bi and p+ = PB2 + PBi ■ The 
outcomes of the measurement {x-,p+} are classically 
stored by Charlie. 

Note that, after the Bell measurement, the conditional 
state of the remaining modes Ri, R2,Ci and C2 is Gaus- 
sian with CM of the form Il3ll 



VR1R2 

XT 



X 



(4) 



where Vr^r^ describes Ahce-Bob remote modes, VdCa 
describes the certification modes in Charlie's hands, and 
X accounts for cross-correlations. From the blocks, 
VrjRj and VC1C2: one computes the partially-transposed 
matrices, Vo^'p and Vf^'r , from which one can derive 
the minimum symplectic eigenvalues [J, |13| 
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The certifying condition of Eq. ^ and the definition of 
logarithmic negativity E^^^ = max{0, — In 277^^5^} [l^, 
17[ imply that 



R1R2 

N 



>E' 



C1C2 

N 



> 0. 



(6) 



In other words, performing the entanglement swapping 
protocol with a certifying state allows to entangle both 
the C1C2 and the R1R2 bipartite subsystems, with the 
latter one more entangled than the former one. There- 
fore the detection by Charlie of entanglement in the cer- 
tification subsystem C1C2 guarantees the generation of 
remote entanglement, i.e., between the remote modes, Ri 
and i?2, in Alice's and Bob's stations. This is the key- 
idea of the protocol. Charlie can locally test the presence 
of entanglement between the certification modes. If this 
is detected, then Charlie is sure that potential entangle- 
ment is also present between Alice and Bob. This remote 
entanglement will be definitely established once that the 
Bell- measurement results are communicated to 

Alice and Bob. 

Therefore in the third step of the protocol, soon af- 
ter the Bell measurement, Charlie detects the certifica- 
tion modes Ci and C2 , by measuring the quadratures or 
other entanglement witness operators. After the success- 
ful detection of C1C2 entanglement, Charlie communi- 
cates the Bell-measurement outcomes to Alice 
and Bob. As a result, the two users now share two en- 
tangled modes, Ri and R2, which are described by a 
Gaussian state with CM Vr^Rj and non-zero first mo- 
ments (Irj^ and (Irj . These moments depend on the pair 
{x-,p+} and are known to Alice and Bob. Thus, the 
two users may decide to delete cIr^ and dp^ by applying 
displacements [13] . More simply, they store the values of 
the first moments and account for them in the classical 
post-processing of data. The previous steps can be re- 
peated many times, with Alice and Bob collecting a set 



of entangled modes {r\^\r^^}, {R^' ^R'^'"'}, ■ ■ • , where 
each pair {r'^\R2^} has the same CM VrjR2 but differ- 
ent first moments {dp^\d^'^^}. 

Application to optomechanical systems. As a first rele- 
vant application of the proposed protocol, we consider the 
case where the initial tripartite systems are two identical 
cavity optomechanical systems [see Fig. ([TJd)]. 

Each system is formed by a micro-mechanical resonator 
playing the role of the remote mode, coupled to two inde- 
pendently driven optical cavity modes with different fre- 
quencies [1^, whose output field are sent to Charlie's site 
(these corresponds to the Bell and certification modes). 
In such an example, the proposed protocol allows to gen- 
erate an entangled state between two remote mechanical 
resonators which can be used for tests of quantum me- 
chanics on massive objects at the nanogram level, and 
whose entanglement can be certified using only standard 
optical homodyne measurements at Charlie's site. In 
Fig. ([TJ;) we have schematically described the optome- 
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chanical system as a bi-chromatically driven Fabry-Perot 
cavity with a movable micromirror, but the present treat- 
ment applies to a generic system in which the optome- 
chanical coupling is linear in the resonator position (e.g., 
see Refs. [HHa). 

The Hamiltonian of the optomechanical system is 



H = He 



M 



OM 



IN, 



(7) 



where Hq — ^ hujkal.ak with fc = b, c is the Hamiltonian 
of the two driven cavity fields, and Hm = fiw,n(p^ + <f)/2 
is the Hamiltonian of the mechanical resonator with mass 
m, resonance frequency Wm, and dimensionless quadra- 
tures p and q with [q,p\ — i. The third term represents 
the optomechanical interaction Hqm = ^hq'^Gokaj.ak, 
with coupling constants Gok — {i^Jk/ L)^/hpmijJm (for 
a Fabry-Perot cavity of length L with a movable mi- 
cromirror). The last term describes cavity driving, 
Hm ~ ifiJ2£kol.exp{—iiJi^^kt) + c.c, with driving rate 
|£fe| = y^2KkPk/lT^'L,k, where cjL,fc is the frequency of the 
input laser, Pk is its power, and Kk is the cavity mode loss 
rate through its input port. The above Hamiltonian is 
valid when the scattering of photons of the driven modes 
into other cavity modes and also between the two chosen 
modes is negligible: this is valid when Wm is much smaller 
then the free spectral range of the cavity 0] . 

The full dynamics of the cavity optomechanical system 
involves fluctuation-dissipation processes, which can be 
taken into account by adopting quantum Langevin equa- 
tions (QLEs) [2^. If both cavity modes are intensely 
driven and the system is in its stability region, the latter 
reaches a semiclassical steady state with large intracav- 
ity amplitudes for both cavity modes, Q;s,fc, and a shifted 
equilibrium position of the resonator, qs. The fluctua- 
tions of the various operators around such a steady state 
obey a quantum dynamics well described by a lineariza- 
tion of the QLE. This corresponds to writing an oper- 
ator O as the sum of its semiclassical value plus quan- 
tum fluctuations O ^ Os + SO ^. The linearized QLEs 
have a compact form in terms of the quadrature fluc- 
tuations of the mechanical mode {Sq, dp} and the two 
cavity modes, 6xk and Sijk, defined by the relation 5ak = 
{Sxk+iSijk)/ V^- In fact, the QLEs can be written as u = 
Au + h, where u = (Sq,5p,Sxh,Syh,5xc,6yc)~'^ . In this 



formula, fi = (0, v2Ki™, v2Ky\ 
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a noise vector, where ^ is the thermal noise affecting the 
resonator, and {x™, y^} are optical input noises. Finally, 
A is the drift matrix of the system quantum fluctuations, 
which is given by 
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where 7™ is the mechanical damping rate, Gk = 
^/2GokCts,k are the effective optomechanical coupling be- 
tween the fluctuations, and — LOk — uj-^^k — Gofers are 
the effective detunings of the two cavity modes. 

We are however interested in the traveling output fields 
associated with the cavity modes. These fields are sent 
to Charlie for the Bell and subsequent certifying mea- 
surements. Considering the cavity output is advanta- 
geous also for achieving larger and more robust entan- 
glement: in fact, by properly choosing the central fre- 
quency and the bandwidth of the output modes, one can 
optimally filter the frequencies carrying the proper op- 
tomechanical quantum correlations [y, |21[ . The selected 
output modes are defined by the following bosonic anni- 
hilation operators a|' 



Ri Ro 



^^\t) = // hkit - s)a^."'(s)ds, where 



= v2K(5afe(i) — a)J^(t) is the standard input-output 
relation for the optical field, and hkit) is a causal filter 



function defining the output modes [21|. In fact, a^^ is a 



standard photon annihilation operator, implying the nor- 
malization condition J \hk{t)\'^dt = f. A simple choice 
is taking hk{t) = ^j2/TkQ{t) exp[— (I/rfc + iVLk)t] where 
Q{t) is the Heaviside step function, f /r^ is the bandwidth 
of the filter, and VLk is the central frequency (measured 
with respect to the frequency of the corresponding driv- 
ing field). 

We are interested in the stationary state of the quadra- 
ture fluctuations of the mechanical resonator and the two 
selected output modes. This is a zero-mean Gaussian 
state since the system is driven by Gaussian noises and 
the dynamics is linearized around a semiclassical point. 
This tripartite Gaussian state, for both Alice's and Bob's 
sites, represents the input state of our entanglement 
swapping protocol. It is fully characterized by the CM 
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ij 
out 



°"*(cx)))/2, where 
whose form can be 



(u°"'(oo)m°^'(oo) + M°"*(oo)^t' 
= {5q, Sp, xf, yl'\ xf, yfV, 
explicitly evaluated using the formalism of Refs. |d l2fj. 

The system parameters must be chosen in such a way 
that the tripartite system formed by the mechanical res- 
onator and the two cavity output fields is in a station- 
ary state satisfying the certifying condition of Eq. ^ 
which, as noticed above, implies the presence of nonzero 
entanglement between the mechanical resonator and the 
output mode employed for the Bell measurement, and 
also between the two optical output fields. The certify- 
ing condition is satisfied if we appropriately choose the 
detunings and filter the appropriate output modes. In 
fact, the optomechanical entanglement is largest when 
we drive the cavity Bell mode with a blue-detuned laser 
(Ab < 0) and the cavity certifying mode by a red-detuned 
laser (Ac > 0). The behavior of the system in an experi- 
mentally feasible parameter region is shown in Fig. [D^a), 
where the resulting log-negativity En between the two 
remote mechanical resonators is plotted as a function 
of the normalized inverse bandwidth of the Bell output 
mode TbUJm, and of the cavity mode bandwidth n/LUm 



(equal for the two modes). Figs. [l{b),(c) show E^^ 
of the remote mechanical resonators (blue circles) and 



C1C2 



of the certification modes (red plus) along lines 1 
and 2 of Fig.[2][a), respectively. In Figs. [2jb),(c), the pa- 
rameter regions where the condition i?^^*^^ > E'^'~^ > 
of Eq. © is satisfied correspond to the situation where 
the initial optomechanical tripartite state satisfies the 
certifying condition. 
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FIG. 2: (Color online) (a) Contour plot of the logarithmic 
negativity En of the remote mechanical resonators with re- 
spect to the normalized inverse filter bandwidth Ttujm, and the 
normalized cavity decay rate K/a;,„. (b) and (c) show En of 
the mechanical resonators (blue circle) and of the certification 
modes (red plus) along line 1 and line 2, respectively. Param- 
eters are: Cavity length L = 1 mm, mechanical resonator 
with frequency ujm/^-K = 10 MHz, quality factor Qm = lO'', 
and mass m — 10 ng, giving coupling constants Got ~ 1 KHz. 
Optical modes have detuning values Ac — — Ab — cJm, while 
the output modes have central frequencies f^b — — i^m and 
Qc ~ ij^m, the relation between the two inverse bandwidths 
being Tc = Tb/5. Input laser powers are Ph = 4.5 mW and 
Pc ~ 5.0 mW, with Xh ~ Xc ~ 810 nm. Reservoir tempera- 
ture is r = 0.1 K. 



In conclusion, we have proposed an extension of en- 
tanglement swapping protocol where the generation of 
entanglement between two remote systems can be certi- 
fied by local measurements performed soon after the Bell 
measurement. We have then shown that this protocol 
can be applied for entangling two distant mechanical res- 
onators, which are part of two optomechanical systems. 
The resulting mechanical entanglement can be success- 
fully certified by applying local optical measurements. 
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of the two subsystems 
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Supplementary Material 
Conditional correlation matrix 

The system is initially composed of two distant, in- 
dependent, tripartite continuous variable (CV) systems, 
RiBiCi in Alice's hands, and R2B2C2 in Bob's hands. 
As a consequence, the total initial state is given by ten- 
sor product of the states prepared by Ahce and Bob, 
Pin = Pi ® P2, which in the Wigner function formalism, 
can be expressed as the product of the Wigner functions 



(9) 



where Offe, /3fc, and 7^, with fc = 1,2, are the complex 
amplitudes related to the phase space variables, e.g. we 
have Uk = (xfc -I- ipk)/V2. 

Then, both Alice and Bob send two of their modes, the 
Bell modes Sfc, and the certifying modes Cfc, to Charlie, 
while holding the remote modes Rk in their hands. At 
the next step, Charlie performs the CV version of the 
Bell measurement on the two Bell modes Bi and B2: he 
mixes Si and B2 with a balanced beam splitter giving 
the following output modes 



x± 



XB2 ± ^Bi 



P± 



PB2 ±PBi 



(10) 



V2 ' V2 ' 

which are then both subjected to homodyne detection 
so that the quadratures X- and are measured, with 
outcomes i_ and p+ and the associated probability 
P{x-,p+). Such a measurement collapses the system 
from a hexapartite to a fully correlated quadripartite 
system in a state which is conditioned to the measure- 
ment outcomes. In the Wigner function formalism this 
is expressed by an integration over a complex amplitude 
mixed by the beam splitter, weighted by the measure- 
ment results probability, i.e., 

W^con = / dWi(ai,7i,/3)W^2(a2,72,/3* -^*) , 

(11) 

where /? = ipj^ ~ x- is the complex variable represent- 
ing the measurement outcomes, and (3 is the unobserved 
variable which is integrated out. The variable (3 is related 
to the Bell modes variables mixed by the beam splitter, 
/3± = (P2 ± /3i)/\/2, with (i± = (a;± + ip±)/V2, by the 
following relation /3 = [x+ -\- ip+ — {x^ + ip-)]/2. 

An equivalent description of the conditional state 
is provided by its symmetrically ordered characteris- 
tic function, which is the complex Fourier transform of 
Wigner function $(A) = J"[M^(q!)] (see e.g. Ref. [jy]). Us- 
ing Eq. pT|) one gets 

*co„(Ai,A2,Aii,/i2) (12) 

= ~7^7ST / d^f'2$l(Al,Ail,l'2)'J'2(A2,M2,i^2) , 

where Ai, A2, /ii, /^2, ^^i, and V2 are the complex variables 
conjugated to the Wigner function variables ai, a2, 71, 
72, /3i, and /32 respectively. 

In the case considered in the Letter, the initial tri- 
partite CV states at Alice and Bob sites are Gaussian, 
i.e., with symmetrically ordered characteristic function 
$j(k) = exp{-k'^Vjk/2 -f idjk}, j = 1,2, with covari- 
ance matrix (CM) 



DJ B, E, 
fJ EJ C, 



(13) 
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(the blocks Hj , Bj , Cj , Dj , Ej , and Fj are 2x2 subma- 
trices), vector of the real- valued system variables k, and 
displacement vector of the state d^- , which can be taken 
equal to zero (displacements do not affect entanglement). 
Inserting these expressions into Eq. (|12l) and performing 
the Gaussian integrals, one can express the effect of the 
Bell measurement performed by Charlie as a relation be- 
tween the CM of the input and output Gaussian states. 
In fact, the initial state in possession of Alice and Bob is 
fully specified by the input 12 x 12 CM 



(14) 



with Vj, j = 1,2, specified in Eq. After the Bell 

measurement, the quadripartite system R1R2C1C2 is in 
a conditional Gaussian state with 8x8 CM 



V 



R1R2 

V, 



C1C2 



(15) 



where the CMs VrjRj and VciC2 of the two bipartite 
subsystems R1R2 and C1C2, and the matrix expressing 
their correlations X are given by 



X 



Ri DiZM iZD7 
D2M iZD| 

EiZM-iZE]" 
E2M-IZE7 



DiZM ^Dj 



R2 



Ci 



D2M 



EiZM-^E 



^2 - E2M-IEJ 

Fi - DiZM-iZEi D1ZM-IE2 
E1ZM-ID2 F2-E2M-ID2 



(16) 



where we have introduced M = Bi + ZB2Z and Z = 
diag[l, —1]. The reduced state of the remote site system 
R1R2 possesses a nonzero displacement vector given by 



2D[ZM^O 



dR, = 2D^M^O 



(19) 
(20) 



where O = (a;_,p_ 
outcomes. 



Y is vector of the Bell measurement 



Symplectic eigenvalues 



We can get simplified expressions by adopting a stan- 
dard form for the CMs of the initial tripartite Gaussian 
states at Alice and Bob sites, Vi and V2. In fact, by 
assuming non-singular 2x2 blocks, the two CMs can al- 
ways be simplified to the standard form of Ref. ^ via 
local unitary operations at Alice, Bob and Charlies sites 
(therefore not affecting the initial entanglement). This 
standard form is given by 



r 




d 


I 


/' 




r 


d' 


/" 


/"' 


d 




b 


e 






d' 


b 




e' 


f 


r 


e 


c 






/"' 


e' 




c 



V = 



This means that we can put the 2x2 blocks of Eq. ([13)) 
in the form ~ r^,!, B^: = bki, = Cfel (I is the 2x2 



(21) 



identity matrix), — diag[c?A;, d'j,], E^ = diag[efe,( 



and 



Ffc 



fk f'k 
fii fill 
Jk Jk 



where fc = 1, 2. As a consequence, the CM of the reduced 
states of R1R2 and C1C2 of Eqs. (I16p and ([TT)) assume 
the block-diagonal form 



VriR2 — 



V, 



C1C2 



bi+b2 
did2 



fcl+b2 



Cl 



bi+b2 



eie2 
bi+b2 



Cl 



bi+b2 



bi+b2 



61+62 



61+62 



did2 
61+62 



r2 



61+62 



61+62 



C2 



61+62 



f2 



^1^2 
' 61+62 



61+62 



61+62 



C2 



(22) 



61+62 



(23) 



One then computes the partially-transposed matrices, mum symplectic eigenvalues r/p^^ and 77(-^^_^ which allow 
^RiR2 ^CiC2' from which one can derive the mini- 
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to quantiiy the remote and certifying entanglement. In 
the general case of different tripartite states at Alice and 
Bob sites the explicit expressions of r]^ ^ and ry^ are 
cumbersome, but they considerably simplifies in the sym- 
metric case of identical initial resources at the two sites, 
Vi = V2, considered in the Letter. In fact, one gets 
^RiR2 ~ b~^^/de^VRB and ij^^^^ = b~^^/det Vbc, where 
Vrb and Vbc are the CMs of the reduced Gaussian in- 
put state of the bipartite subsystems RiBi {R2B2), and 
BiCi {B2C2), respectively. Using the expressions for the 
= [4VdctVRBj^ ^ 



put in bipartite standard form. 



purities, /i^ 

and /Xg = 
Letter, 



2Vdct B 



Aiec = [4VdetVBcJ , 
one arrives at Eq. (5) of the 



V-RB = 



Vbc - 



r d 

r d' 
d b 

d' b 

b e 

b e' 
e c 



(25) 



(26) 



Therefore one can apply the derivation of Ref. [3| to both 
the RB and BC subsystems. 



^'RiR2 ^ 



2Mrb 



(24) 



This derivation of the partial transpose minimum sym- 
plectic eigenvalues practically coincides with that of 
Ref. Q, which refers to entanglement swapping start- 
ing from identical bipartite Gaussian states in standard 
form. In fact, the tripartite standard form of Ref. 2] im- 
plies that the BC and RB subsystems are simultaneously 
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